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ABSTRACT 

Anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) are recognized as the most 
promising magnetar candidates, as indicated by their energetic bursts and rapid spin-downs. It is 
expected that the strong magnetic field leaves distinctive imprints on the emergent radiation both by 
affecting the radiative processes in atmospheres of magnetars and by scattering in the upper magneto¬ 
spheres. We construct a self-consistent physical model that incorporates emission from the magnetar 
surface and its reprocessing in the three-dimensional (3D) twisted magnetosphere using a Monte Carlo 
technique. The synthetic spectra are characterized by four parameters: surface temperature kT, sur¬ 
face magnetic field strength B , magnetospheric twist angle A tf>, and the normalized electron velocity 
/3. We also create a tabular model (STEMS3D) and apply it to a large sample of XMM-Newton spectra 
of magnetars. The model successfully fits nearly all spectra, and the obtained magnetic field for 7 out 
of the 11 sources are consistent with the values inferred from the spin-down rates. We conclude that 
the continuum-fitting by our model is a robust method to measure the magnetic field strength and 
magnetospheric configuration of AXPs and SGRs. Investigating the multiple observations of variable 
sources, we also study the mechanism of their spectral evolution. Our results suggest that the magne¬ 
tospheres in these sources are highly twisted (Ac/) > 1), and the behavior of magnetospheric twisting 
and untwisting is revealed in the 2002 outburst of IE 2259+586. 

Subject headings: radiation mechanisms: nonthermal — stars: magnetic fields — stars: neutron - 
X-rays: stars 


1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) and soft gamma- 
ray repeaters (SGRs) form an intriguing class of iso¬ 
lated neutron stars (NSs). They share the similar ob¬ 
servational properties: (1) relatively slow spin peri¬ 
ods (P ~ 2 — 12 s); (2) large spin-down rates (P ~ 
10 -14 — 10 -10 s s _1 ); (3) high persistent X-ray lumi¬ 
nosity (Lx ~ 10 34 — 10 36 erg s _1 ), which significantly 
exceeds their rotation power for most sources; and (4) 
different types of X-ray variability; i.e., long-term X- 
ray flux enhance men ts and short energeti c bursts (see 
lMereghettill2008L and lRea fe Esuositoll201 ll for recent re¬ 
views). These behaviors can be understood within the 
magnetar scenario - extremely strong magnetic fields 
of some N Ss givin g rise to the observed exo tic man¬ 
ifestations (Thompson fc Duncanl ll995l . 199 6). There 
are now 28 magnetars or magnetar candi dates known Q 
(jOlausen fe Kaspill20l4 IZhou et aI1l2014fl . 

Although recent multi-wavelength observations of 
magnetar activities significantly widened our views, the 
most essential knowledge on magnetars is provided by 
the emission properties of their X-ray radiation. The X- 
ray spectra of SGRs and AXPs below 10 keV are soft and 
can be empirically described by a power-law (PL) model 
with photon index r ~ 2 — 4 and a ther mal com ponent 
with temperature kT ^0.5 keV ijMereghetti et al.ll2005l 
2006). Also note that it has been revealed that some 
magnetars ex hibit a very hard spectrum (r ~ 1) above 10 
keV (see, e.g. JKuiper et al .112004[2006tIclen Hartog et all 

1 http://www.physics.mcgill.ca/$\sim$pulsar/magnetar/main. 


I2008allbh . iBeloborodovI (20131 1 suggested that the hard 
X-ray emissions attribute to surface emissions scattered 
by the relativistic outflows around magnetars, which also 
make a non-negligible contribution be low 1 0 k eV. Study¬ 
ing th e spectra of persistent emission, iMarsden fe Whitel 
(1200 lfl reported that the hardness of the spectra increases 
with the spin-down rate v . With more accumulated data, 
iKasni fc Bovdstuiil ( 2010 1 showed that the spectral in¬ 
dices in quiescence are not only correlated with v but 
also the magnetic field B inferred from timing analy¬ 
ses. Magnetars occasionally emit short (~ 0.1 s) bursts, 
during which the plausible line fe atures ha ve been also 
detec ted in some sources (e.g., iStrohmaver fc Ibrahiml 
l200d) . If such features are due to the proton cyclotron 
absorption, the inferred magnetic field somehow deviates 
from the value estimated from the spin-down rate (e.g ., 
iWoods et al.l 120051 : iTiengo et all 120131 : 1 An et al.1 12014 1. 
An independent method is thus required to evaluate the 
magnetic field of magnetars. 

In the magnetar framework, the X-ray radiation 
can be physically interpreted as the NS surface emis¬ 
sion scatte red by the plasma in the magnetosphere 
(see, e.g.. iThompson et al.1 120021 iLvutikov fe Gavriill 
2006|: |Fernandez fe T hompsonl 120071: iTong et al.l I2010H . 
Thompson et al.1 (l2002ll suggested that the internal mag¬ 
netic field of magnetars is highly twisted and anchored 
in the highly conducting crust of an NS. In starquakes, 
the strong internal magnetic field deforms the crust, 
and as a result the footpoints of the external magnetic 
field lines are displaced, i.e., the magnetic field out¬ 
side the star - magnetosphere is also twisted. Alterna- 
iitmi'ively, a twisted magnetosphere can also be aroused by 
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continued footpoint motions (|Beloborodovl l2009h . The 
twisted external magnetic field can support electric cur¬ 
rents that are much larger th a n the Goldreich-Julian cur - 
rent dThompson et alJ l2002t iLvutikov fc Gavriil l200 6j. 
The charge carriers provide large optical depth to res¬ 
onant cyclotron scattering (RCS) so that soft photons 
can be Compton up-scattered and form a PL high- 
energy tail. Both one-dimensional (ID) and three- 
dime nsional (3D) RCS models had been developed (see 
e.g. , ILvutikov fc Gavriil 120061 : iFernandez fc Thompson! 
2007), and offer ed promising fit s for the X-ray sp ectra 
of SGRs/AXPs (|Rea et alJl2008t iZane et all 1200911 . Be¬ 
cause the non-resonant scattering contributes a negligi¬ 
ble fraction, both ID and 3D RCS models only take the 
resonant scattering into consideration. For ~ keV pho¬ 
tons, the resonant scattering occurs in the layer where 
the magnetic field decays to ~ 10 11 G regardless of the 
NS surface magnetic field. As a consequence, the es¬ 
sential information of SGRs/AXPs - the strength of the 
magnetic field —is not incorporated in both ID and 3D 
RCS models since the scattering process is insensitive to 
the magnetic field below its cyclotron frequency (> 100 
keV for a magnetar magnetic field). 

Besides the magnetospheric scattering, the surface 
emissions are also distorted in strongly magnetized (B > 
10 13 G) NS atmospheres due to the effect of vacuum 
polarization and proton cyclotron resonances, by af¬ 
fecting the interactions between the photons and the 
plasma (e.g., IHo fc bail 12003 [CDdl [2(lo3 . Thus, the 
spectral profile of surface emission strongly depends 
on the magnetic field. Taking into account the com¬ 
bined effect s of the mag netar atmosphere and its mag¬ 
netosphere, iGiiver et al.l (|2007l ) developed the ID sur¬ 
face thermal emission and magnetospheric scattering 
(STEMS) model, which can successfully fit the soft 
X-ray spectra of both steady and variable magnetar s 
(IGiiver et a,l.ll20 07. 2012: INg et al 11201 IlibbTet, al.ll2012H . 
Additionally, the STEMS model uncovers the informa¬ 
tion concering the magnetic field. However, because the 
STEMS model treats the scattering region as a plane- 
parallel slab, it cannot offer information about the ge¬ 
ometry of the magnetosphere. 

The long-term monitoring observations revealed that 
the X-ray emi ssions of most magnetars are variable 


•ay 

le.g-. lRiallMl . in the extreme cases, namely the tran¬ 
sient magnetars (e.g., XTE J1810-197, and IE 1547.0- 
5408), the overall flux can be enhanced up to three 
orders of magnitude brighter than their quiescent 
level (e.g.. Ilbrahim et alJ 12004 iMereghetti et al.l 120091 : 
iKaneko et al.l I20l7ib Studying X-ray behaviors in dif¬ 
ferent flux levels would help us to explore the na¬ 
ture of magne tar, i.e., the evolution of magnetospheres 
(lBelobor odovll20 (Ml and the magneto-thermal e volution 
of NSs (jPons fc Rea I l2012t iVigano et al.l 201 3fl . Based 
on XMM-Newton observations, IZhn et al.N2008h found 
that the afterglow of the IE 2259+586’s 2002 outburst 
followed by a PL decay and the hardness is strongly cor¬ 
related with the flux. Alternatively, the flux relaxation 
of XTE J1810-197’s 2003 outburst is best described by 
an exponential decay and t he surface temperatu re be¬ 
came cold in the meantime (jGotthelf fc Halp'ernl l2005il . 
Diver se observational phenomena (see iRea fc Espositol 
l20TH and reference therein) suggest that X-ray activi¬ 


ties could be driven by d ifferent mechanisms, ei ther the 
twisted magnetosphe res dThompson et al.l I2002D or the 
deep crustal heating (jLvubarskv et al .Il2002f~ 

In this work, we carry out 3D Monte Carlo simulations 
of the emitted photons from the surface of a highly mag¬ 
netized NS propagation in a twisted magnetosphere. The 
physical bases of our model and the Monte Carlo method 
are described in Section 2. The model properties are pre¬ 
sented in Section 3. We also numerically calculate model 
spectra and create a tabular model, which is further ap¬ 
plied to the X-ray spectra of magnetars (Section 4). The 
ultimate goal of this paper is to better understand the 
surface and magnetospheric properties of magnetars and 
their evolutions using our model. In Section 5, we discuss 
the spectral modeling results and their implications. 

2. MODELS 

In this section, we first outline the physical processes 
in a strongly magnetized atmosphere, and then scatter 
the surface emission in the 3D twisted magnetosphere. 
The gravitational redshift effect is also considered in our 
model. 


2.1. Surface emission of magnetars 

In strong magnetic fields (B > 10 13 G), the vacuum po¬ 
larization, a quantum electrodynamics phenomenon, can 
affect radiative transfer in magnetized plasmas by mod¬ 
ifying polarization modes and the opacities of the nor¬ 
mal m odes (e.g., IZane et al.l f2001fc IHo fc Laill2003 lOzell 
2003). The proton cyclotron resonance is another phe¬ 
nomenon that affects the propagation of photons in a 
magnetized plasma. Following the methods discussed in 
lOzell (1200 ll . 12001 . the radiative equilibrium models are 
constructed for the fully ionized hydrogen plasmas, tak¬ 
ing into account the effects of vacuum polarization and 
ion cyclotron lines. We use a modified Feautrier method 
for the solution of the angle- and polarization-mode de¬ 
pendent radiative transfer problem and ensure radiative 
equilibrium with a temperature correction scheme based 
on the Lucy-Unsold algorithm. The surface emission 
spectrum can be defined by just two parameters: the 
effective temperature of atmosphere kT and the surface 
magnetic field strength B. 

2.2. Resonant cyclotron scattering in 3D twisted 
magnetosphere 

The Monte Carlo techniques are quite suitable to han¬ 
dle photon scattering in complicated 3D configurations 
and the 3D RCS Monte Carlo code has already been 
built b y [Fernandez fc Thompson! (120071 ) and lNobiii et alJ 
(|2008l) . Both of these works considered axisymmetric 
globally twisted magnetospheres and seed photons as be¬ 
ing the canonical blackbody (BB) emissions. Neverthe¬ 
less, these codes allow for an arbitrary distribution of 
seed photons, velocity distribution of the charged parti¬ 
cles, and magnetic field geometry. 

We build our Monte Carlo code fo llowing an ap proach 
that is similar to the one discussed in lNobili et al.l (120081) . 
The fundamental difference is the description of seed pho¬ 
tons. (1) The original BB emissions are distorted af¬ 
ter travelling through strongly magnetized atmospheres, 
most ly carried by extraordinary mode photons (l Ozell 
120031 ). In our work, this surface emission is adopted as 





























































































STEMS3D 


3 


the seed photons, which are homogeneously injected from 
the stellar surface with a wavevector pointing in the ra¬ 
dial direction. Note that under the influence of an ultra¬ 
strong magnetic field, the surface emissions are expected 
to be anisotr opic and their temperature varies with 
latitu de (e.g., iHevl fc Hernmiistl 1 1 9981 iBernardini et, al.1 
unnD. However, the study on the anisotropic surface 
emissions is beyond the scope of this paper and will 
be reported elsewhere. (2) The globally twisted mag¬ 
netosphere has a self-similar and force-free construction 
(j x B = 0), which is uniquely characterized by the net 
twist angle of the field lines anchored close to the two 
magnetic poles, A<j>. (3) The current density j is deduced 
from the Ampere’s law (V x B = (47r/c)j) un der the 
f orce -free assumption. Following the work in lNobili et al.1 
in, we consider the effects of bulk velocity and also 
thermal velocity distribution of charge carriers, assuming 
a ID Maxwellian distribution at a given temperature su¬ 
perimposed to a bulk motion. In order to reduce parame¬ 
ters in the model, the electron temperature kT e is derived 
from the bulk velocity (and further halved) by assuming 
the equipartition between thermal and bulk kinetic en¬ 
ergy. As a result, the number density of electrons (n e ) 
is a function of the twist angle and the average electron 
velocity. (4) We neglect both the non-resonant scatter¬ 
ing and electron recoil effects, which are not important 
in the soft X-ray band. The spectrum is calculated in 
the non-relativistic regime with t he s impli fied r esonant 
cross-section (Equation (10) in lNobili et al.il20(M) . which 
is independent of magnetic field and frequency. 

As photons propagate in the 3D twisted mag¬ 
netosphere, we check whether they can freely 
escape out of the m a gnetos p here in eve r y ste p 
(Fernandez fc Thompsonl 120071 iNobili et al.l 12008 ). 
If the photon frequency together with moving direction 
still allow resonant scattering, as soon as the integrated 
scattering depth r s satisfies r s > —InU, where U is 
a random number in the range [0,1], a scattering is 
triggered, and the photon frequency, direction, and 
polarization are updated. Then, the photon continues 
to travel in the magnetosphere until another scattering 
occurs o r it escapes to infinity (observers). We also 
refer to INobili et al.l (l2008h for more details on the 
scattering process. The spectra are calculated by 
evolving 5,000,000 photons and characterized by four 
parameters: surface temperatures kT, magnetic field 
strength at the poles B, magnetospheric twist angles 
A</>, and the normalized electron velocity (3 = v/c, where 
c is the speed of light. An example spectrum (kT = 0.3 
keV, B = 10 14 G, A(f> = 1.0, and /3 = 0.3) is shown in 
Figure Q] 

2.3. Gravitational redshift effect 

Because of the strong gravitational field of NSs, the 
gravitational redshift effect is important and should be 
corrected for the photons close to NSs, but is negligible 
for the photons that are far away from NSs (> several 
times of the NSs radii). According to General Relativity, 
the photons experience a redshift in gravitational fields 
by a factor of 1 + z = , 1 2QM assuming Schwarzschild 

V 1 - Rc 2 

geometry. In principle, the gravitational redshift should 
be corrected for the surface emission, before the magne¬ 
tospheric scattering. Since the scattering optical depth 
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Figure 1. The BB (kT = 0.3 keV, blue solid line) and the mag- 
netar surface (MS, kT = 0.3 keV and B = 10 14 G, red dotted line) 
emissions are scattered in a 3D twisted magnetosphere (A <j> = 1.0 
and 0 = 0.3), resulting in the spectra of the 3D RCS (blue dashed 
line) and the STEMS3D (red dotted-dashed line), respectively. 
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Figure 2. Example spectrum (kT = 0.3 keV, B = 10 14 G, A <j> = 
1.0, and 0 = 0.3) with the gravitational redshift correction (z = 
0.306) performed before (blue plus) and after (red diamond) the 
magnetospheric scattering. 

is inde pendent of frequency in a self-similar magneto- 
sphere dThompson et all 1200 2; iFernandez & Thompsonl 
I2007D . it allows us to correct the spectra after the mag¬ 
netospheric scattering. In order to test this hypothesis, 
we calculated the spectra for kT = 0.3 keV, B = 10 14 
G, A(f> = 1.0, and /3 = 0.3 with the gravitational red¬ 
shift correction of z = 0.306 (corresponding to an NS 
with mass 1.4 M 0 and Rns = 10 km) performed before 
(red line in Figure [2]) and after (blue line) the magneto¬ 
spheric scattering, and find that both spectra are identi¬ 
cal. Therefore, it becomes feasible and convenient to do 
the correction in X'SPEC (after scattering, §4). 

3. MODEL PROPERTIES 

The emerging radiation from the surface of a highly 
magnetized NS is expected to be modified by both vac¬ 
uum polarization and proton cyclotron resonance. At 
magnetar-type field strengths, the resulting spectrum is 
spectrally harder than a Planckian and would exhibit 
proton line features in the soft-X-ray band (Figure [TJ. 
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Figure 3. Angle-averaged optical depth (r) vs. A<p and different 
values of /3: 0.2 (solid), 0.3 (dotted), 0.5 (dashed), and 0.8 (dashed- 
dotted). 

The emerging photons further gain energy by multi¬ 
ple scattering in the twisted magnetosphere, manifest¬ 
ing themselves as the high energy tail. To illustrate this 
point, we performed the scattering of a BB ( kT = 0.3 
keV) as well as the magnetar surface emission (kT = 0.3 
keV and B = 10 14 G) in a 3D twisted magnetosphere 
(A (j> =1.0 and (3 = 0.3). We present these spectra in 
Figure[T] We also calculated the spectrum of the 3D RCS 
in the same way as the spectrum of STEMS3D model, 
but with a single temperature BB adopted for the seed 
photons. 

Unlike the STEMS model, in which the optica l dept h is 
indep e ndent of the partic le velocity (ILvutikov fe Gavrul 
2006{ iGiiver et al.l 1200711 . the optical depth in the 
STEMS3D is determined by both the degree of magne- 
tospheric twist and the velocity of charged particles. In 
Figure [3] we show that the variation of angle-averaged 
mean number of scattering per photon (r) with respect 
to A</> and (3. Figure [4] illustrates the effects on the spec¬ 
tral shape of varying twist angle and electron velocity. 
It is clearly seen that the level of upscattering of emer¬ 
gent thermal photons increases with A 4> and in particular 
with (3. 

ILvutikov fe Gavrul (12006H pointed out that the deep 
line features present in the surface spectra can be 
smoothed out with a large scattering optical depth. How¬ 
ever, investigating the spectra obtained with a narrow 
velocity distribution (/3 = 0 .5, A0 = 0.1) and large twist 
(A (j> = 1.01. [Fernandez fe Thompsonl ([2007D argued that 
the line strength can be reduced by a factor of ~ 2; how¬ 
ever, a significant line feature remained as a result of an 
optical depth in the 3D twisted magnetosphere that was 
not very large. As already mentioned, the optical depth 
in 3D twisted magnetosphere depends not only on A (j> 
but increases sharply with decreasing (3 (Figure [3]). As 
we show in Figure [j we find that the STEMS3D model 
can indeed wash out the line features if the twist angle is 
large (Ac/) > 1) and the electron speed is slow (/3 < 0.3) 
enough. 

4. APPLICATION TO SOFT X-RAY SPECTRA OF 
MAGNETARS 

We generated model X-ray spectra in the 0.01 10.0 
keV range using the numerical model described in §3. 


We used the following parameter space for our numer¬ 
ical grid: the surface temperature kT = 0.1 — 0.6 keV 
(step 0.1 keV), surface magnetic field strength B = 
10 14 — 10 15 G (step 10 14 G), magnetospheric twist an¬ 
gle A cj) = 0.0 — 2.0 radian (step 0.1 radian), and the 
normalized electron velocity [3 = 0.1 — 0.9 (step 0.1). We 
then employed the routine wftbmd, El, to create a tab¬ 
ular model (named STEMS3D.mod) which can be im¬ 
plemented into the standard package for X-ray spectral 
analysis, XSPEC. This tabular model is further applied 
to the observational data of magnetars. 


4.1. Observations and data analysis 


The main objective of our efforts here is to better un¬ 
derstand the surface and magnetospheric properties of 
magnetars using the STEMS3D model. For this investi¬ 
gation, we selected the sources that were bright enough 
(/o.5—iokeV > 10 -12 erg cm -2 s _1 ) to allow statistically 
significant spectral results and were located in relatively 
low interstellar absorption (nH < 5 x 10 22 cm -2 ) regions 
of our Galaxy. When nH > 5 x 10 22 cm -2 , X-ray emis¬ 
sions below 2 keV are strongly absorbed, making the pa¬ 
rameters of the STEMS3D model unconstrained, even for 
the bright source SGR 1806-20. The flux of the serendip- 
itously discovered magnetar, 3XMM J185246.6+003317, 
during its active phase is a few times of 10~ 12 erg cnr 2 
s~ 4 . However, it was only detected with EPIC-MOS 
(|Zhou et al.l 120141) , which has a smaller effective area. 
We do not include these spectra owing to the low signal- 
to-noise ratio (S/N). 

Our sample includes bright persistent magnetars 
that exhibited only subtle long-term flux varia¬ 
tions even when they experienced glitc hes or emit¬ 
ted energetic X-ray bursts 4U 01 42+61 (IGavriil e t al 


2011 ): 1RXS ,1170849. 0-400910 (ISasmaz Mus fc GdfHis 
ili : IE 1841 -045 (ISasmaz Mus et al. I 120141) : and 


SGR 1900+14 (llsrael et al.l I2008D . On the other 
hand, two bright persistent sources: IE 2259+586 and 
IE 1048.1-5937 showed long-term X-ray flux enhance- 
ments in connection with X- r ay b ursts and glitches 
(|Kast>i et al. I 120031 : IDib et ahl 120091) . Over the last 
decade, transient magnetars have emerged as a sub¬ 
group. In quiescence, these sources emit X-rays at levels 
near or below our detection capabilities. Their X-ray 
fluxes are enhanced by a factor of 100 or more at the 
onset of their out burst episodes (IReall20l1) . We selected 
IE 1547.0-5408 (llsrael et all I2008D . Swift J1822.3-1606 


(Sch olz et al.ll2014l) . SGR 0501+456 (Gogus et al 
ILin et al.H20m and XTE J1810-197 rtGiiver et al. 


20H ; 


12007 ) 


to investigate the spectral behavior of transient magne¬ 
tars. Note that the quiescent spectra of transient mag¬ 
netars are excluded even when their fluxes exceed the 
critical flux level. We discuss the implications of these 
spectra and limitations of our model for sources in qui¬ 
escence in §5. We also include in our sample the longest 
XMM-Newton observation of CXOU J171405.7-381031 
in 2010. The observational details of the X-ray data 
used in this paper are listed in Table [I] Note that the 
two short observations of IE 1841-045 in 2002 October 
were performed just two days apart. To improve count¬ 
ing statistics, we combine these two short observations 
to one spectrum. 


2 http://heasarc.gsfc.nasa.gov/docs/heasarc/ofwg/docs/general/\model: 
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Figure 4. Left: model spectra (without gravitational redshift correction) for B = 10 14 G, kT = 0.3 keV, ft = 0.3 and different values of 
A <j>. Right: computed spectra for B = 10 14 G, kT = 0.3 keV, A <ft = 1.0 and different values of ft. 


The data collected with the XMM-Newton EPIC-pn 
instrument are reduced using the Science Analysis Sys¬ 
tem software (SAS) version 12.0.1, and filtered with the 
standard criteria: cleaning for background areas, setting 
FLAG=0 and PATTERN<4. We use the SAS task epat- 
plot to evaluate the pile-up fraction. For observations 
showing severe problems with pile-up (Table [[]), spectra 
are extracted within annulus regions with the inner ra¬ 
dius ~ 10-20" and ARF files are calculated to correct 
the missing part of the point-spread function, thus the 
correct flux level could still be measured through spec¬ 
tral fitting. The spectral response files are created using 
the SAS tasks rmfgen and arfgen. All spectra are re¬ 
binned with the task specgroup to have at least 20 counts 
per bin to enable the use of chi-square statistics and not 
to oversample the instrument energy resolution by more 
than a factor of three. A 2% systematic error is added 
to the data to account for uncertainties in instrumen¬ 
tal calibrations. The spectral analysis is performed in 
the 0.5-7.0 keV energy range for the spectra having a 
low absorption, while in the 1.0-7.0 keV for those being 
highly absorbed. We note that the model-independent 
residuals are detected in the spectra of IE 1547.0-5408 
and CXOU J171405.7-381031 below 1.2 keV; therefore, 
only the spectra in the energy 1.2-7.0 keV are u sed ( Ta- 
ble[l]). All spectra are fitted with XSPEC 12.8.1 (jArnaudl 
119961 ). Because the interstellar hydrogen column density 
(nH) is not expected to vary during different observa¬ 
tions, for each source with more than one pointing, we 
fit all available spectra simultaneously and link the nH 
to have a commo n value, an d the more updated solar 
abundances by iLoddersl lj2003l) are adopted in our work. 

We first attempt to fit all spectra with an absorbed 
STEMS3D model assuming the gravitational redshift pa¬ 
rameter of 0.306. We find that spectra of 4U 0142+61 
and 1RXS J170849.0-400910 clearly require an additional 
hard-X-ray component; therefore, adding a PL compo¬ 
nent significantly improved the fit since the mechanism 
responsible for the hard-X-ray emission likely contributes 
a non-negligible fraction of soft-X-ray emission and nat¬ 
urally affects our fitting (see Figure 0. We fixed the 
photon index T to the values obtained from the fitting 
to the contemporaneous hard-X-ray data. As the spec¬ 
tra from 4U 0142+61, 1RXS J170849.0-400910, IE 1841- 


045, and SGR 1900+14 above 10 keV are relatively sta¬ 
ble and can be fitted with a P L mo del, we adopt a 1' 
of 0.93 (Iden Hartog et al.ll2008aD . 1.13 (Iden Hartog et al. 
2008bf) . 1.32 (|Kuiper et al.ll2006lh and 1.43 (|Enoto et al. 
2010al) . respectively, in addition to STEMS3D in our 
spectral fitting. Because the hard-X-ray component in 
IE 1547 . 0-5408 ex hibited dramatic transient behavior 
(jKuiner et alJlMT^l . we only analyze the XMM-Newton 
observation executed on 2009 February 3, and we used 
a PL index T of 1.41 accordin g to the Suzaku HXD-PIN 
data on 2009 January 28 (|Enoto et al.ll2010bl) . 

4.2. Spectral analysis, spectral evolution, and 
correlations 

4.2.1. Persistent Magnetars 

4U 0142+61 is the brightest magnetar and has 
been one of the most stable X-rav emit ters in X-rays 
(jGonzalez et al.l [20101 : iWang et all I20l4 ). Its longest 
XMM-Newton spectra in the 0.5—7.0 keV band can be 
modeled perfectly with an absorbed STEMS3D plus a 
PL. We find the surface temperature, kT = 0.33 keV, and 
surface magnetic field, B = 5.87xl0 14 G (see Table[2]for 
all other fit details and Figure [6] for the unfolded spec¬ 
trum and best-fit model curve). Note that these surface 
parameters are consisten t with those that were obtained 
using STEMS (iGiiver et alJ I2007D . We find that the 
magnetospheric electrons are non-relativistic (/3=0.21), 
and that the magnetospheric twist angle, A (f> = 1.77, 
that is close to upper end of our parameter space. The 
STEMS3D component contributes more than 98% of the 
total flux below 7 keV. 

Since the first glitch activity detected in 1999 
(|Kaspi et al. I I2000D . sev eral glitches were unveiled i n 
1RXS J170849.0-400910 (ISasmaz Mus fc Gogiis ll20l3l) . 
The source emits persistent radiation in both soft- and 
hard-X-rays. Fitting the 0.5—7.0 keV spectrum, we find 
that the surface temperature is relatively high {kT = 
0.49 keV) and the other parameters are quite close to 
those in 4U 0142+61, i.e., B = 6.03xl0 14 G, A<(> = 1.82, 
and /? = 0.20 (Table0. The STEMS3D component con¬ 
tributes about 88% of the total soft-X-ray flux (Figure 
0 . 

The AXP IE 1841-045 is located at the center of 
the X-ray and radio supernova remnant (SNR) Kes 73 
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Table 1 

Log of the XMM-Newton observations used in this work 


Source Obs No. ObsID Obs Date Mode Net Exposure Energy band 

(ksec) (keV) 


4U 0142+61" 

Obsl 

0206670101 

2004 Mar 01 

Timing 

37 

0.5-7.0 

1RXS J170849.0-400910" 

Obsl 

0148690101 

2003 Aug 28 

SW 

31 

0.5-7.0 

IE 1841-045# 

Obsl 

0013340101 

2002 Oct 05 

LW 

2 

1.0-7.0 



0013340201 

2002 Oct 07 

LW 

4 


SGR 1900+14" 

Obsl 

0305580101 

2005 Sep 20 

FF 

21 

1.0-7.0 


Obs2 

0305580201 

2005 Sep 22 

FF 

20 



Obs3 

0410580101 

2006 Apr 01 

FF 

14 



Obs4 

0506430101 

2008 Apr 08 

FF 

23 


IE 1048.1-5937 

Obsl 

0112780401 

2000 Dec 28 

FF 

4 

0.5-7.0 


Obs2 

014786010U 

2003 Jun 16 

FF 

46 



Obs3 

0164570301 

2004 Jul 08 

Timing 

31 



Obs4 

0307410201 

2005 Jun 16 

SW 

19 



Obs5 

0307410301 

2005 Jun 28 

SW 

16 



Obs6 

0510010601 

2007 Jun 14 

SW 

34 



Obs7 

0654870101 

2011 Aug 06 

FF 

76 



Obs8 

0723330101^ 

2013 Jul 22 

FF 

54 


IE 2259+586 

Obsl 

0057540101 1 " 

2002 Jan 22 

FF 

10 

0.5-7.0 


Obs2 

0038140101 

2002 Jun 11 

SW 

34 



Obs3 

0155350301 

2002 Jun 21 

SW 

18 



Obs4 

0057540201 t 

2002 Jul 09 

FF 

7 



Obs5 

0057540301+ 

2002 Jul 09 

FF 

11 



Obs6 

0203550301 t 

2004 Feb 20 

FF 

4 



Obs7 

0203550601 

2004 Jun 06 

SW 

5 



Obs8 

0203550401 

2004 Jun 22 

SW 

4 



Obs9 

0203550501 

2004 Dec 19 

SW 

4 



ObslO 

0203550701 

2005 Jul 28 

SW 

3 


CXOU J 171405.7-381031 

Obsl 

0606020101 

2010 Mar 17 

FF 

81 

1.2-7.0 

X'l'E J1810-197 

Obsl 

0161360301 

2003 Sep 08 

SW 

7 

0.5-7.0 


Obs2 

0161360501 

2004 Mar 11 

LW 

10 



Obs3 

0164560601 

2004 Sep 18 

LW 

23 



Obs4 

0301270501 

2005 Mar 18 

LW 

36 


SGR 0501+456 

Obsl 

0560191501 

2008 Aug 23 

SW 

34 

0.5-7.0 


Obs2 

0552971101 

2008 Aug 29 

SW 

17 



Obs3 

0552971201 

2008 Aug 31 

SW 

7 



Obs4 

0552971301 

2008 Sep 02 

SW 

14 



Obs5 

0552971401+ 

2008 Sep 30 

LW 

28 


Swift J1822.3-1606 

Obsl 

067228180U 

2011 Sep 23 

LW 

10 

0.5-7.0 


Obs2 

0672282701+ 

2011 Oct 12 

LW 

17 



Obs3 

0672282901 

2012 Apr 05 

LW 

23 



Obs4 

0672283001 

2012 Sep 08 

LW 

20 


IE 1547.0-5408" 

Obsl 

056018110U 

2009 Feb 03 

FF 

49 

1.2-7.0 

SGR 0418+5729 

Obsl 

0610000601 

2009 Aug 12 

SW 

41 

0.5-7.0 

CXOU J 164710.2-455216 

Obsl 

0404340101 

2006 Sep 22 

FF 

40 

1.0-7.0 


Obs2 

0410580601 

2007 Feb 17 

LW 

16 



Obs3 

0505290201 

2007 Aug 19 

LW 

22 



Obs4 

0505290301 

2008 Feb 15 

LW 

19 



Obs5 

0555350101 

2008 Aug 20 

LW 

28 



Obs6 

0679380501 

2011 Sep 27 

LW 

15 



Note. — Mode: operating mode of the EPIC-pn including full-frame (FF), large-window (LW), small-window (SW), and timing 
modes. Net exposure: clean exposure (in unit of kiloseconds) after background flares excluded, energy band: Energy band in which 
spectra are fitted with XSPEC. 

N; Sources’ spectra are fitted by the absorbed (STEMS3D + PL) model with the PL index T fixed (see the text for more details), and the 
others are fitted by a single absorbed STEMS3D model. 

f: Observations suffer problems with pile-up, and spectra are extracted from annular regions centered on the source position. 


(|Vasisht fc Gottlieli 11997ft and has the slowest spin pe¬ 
riod of ~ 11.8 s ( Dib et a,Lll2008l) . The broadband (1-200 
keV) X-ray spectrum can be desc ribed by a BB plus a 
single PL model (|Rea et al.l 120081 b As can be seen in 
Figure [6] the PL component only supplies ~ 8% of the 
total flux below 7 keV; however, it starts to overcome 
the STEMS3D component around 5 keV. The fitting re¬ 
sults indicate that the magnetosphere is modestly twisted 
(A <j> = 0.31) and the electrons have the fastest velocity 
(/3 = 0.28) among all sources studied here, i.e., the small¬ 
est optical depth in the magnetosphere (Figure [3]). The 
model parameters, however, have relative large errors be¬ 


cause of the short effective exposure (~ 6 ks) in addition 
to high interstellar absorption. 

CXOU J171405.7-381031 is associated with the SNR 
CTB 37B and probably with a TeV source HESS J1713- 
381 (lAharonian et al.ll200l . CXOU J171405.7-381031 
was identified as the younge st magneta r with a char¬ 
acteristic age of ~ 950 yr dSato et al.l |2010| ). There 
have been only modest flux variations (less than a fac¬ 
tor of two) reported with two Chandra and one XMM- 
Newton observations du ring 2007 February and 2010 
March (|Sato et all 12010( 1. Currently, no model provides 
a good fit to its X-ray spectra below ~1.2 keV, likely 
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Energy (keV) Energy (keV) 

Figure 5. Spectra of 4U 0142+61 and 1RXS J170849.0-400910 are fitted with the single absorbed STEMS3D model, and the fit residuals 
are shown in panels (a), (b), respectively. Panels (c) and (d) plot the residuals when the same spectra are fitted with the absorbed 
(STEMS3D+PL) model. 


Table 2 

Spectral fit results of 4U 0142+61, 1RXS J170849.0-400910, and IE 1841-045 


Source 

nH 

(10 22 cm" 2 ) 

kT 

(keV) 

B 

(10 14 G) 

A 0 
(rad) 

P 

Flux 

X 2 /dof 


4U 0142+61 

r» '7q - l - U.U3 

U ‘ ' °- 0.05 

p\ qq+LLUcJ 

u.oo 0 03 

r o«+U.t)8 

0,0 - 0.49 

1 77-l-u.uy 
'- 0.06 

0 21 l " u,ui 

11.7 

80.2/122 

1RXS J170849.0-400910 

1 50+°- U5 

0 49+ u - Ui 

u -^_0.02 

6 Q3+ u - ;iU 

O.UO_o 14 

-1 0Q+O.U8 

0.08 

0 20” t " U Ui 

u.^u_o 01 

3.26 

85.2/109 

IE 1841-045 

4 7Q+ U - 81 

' y -1.19 

0 24" i u ‘^ u 

U.Z4_o 07 

^• 1 D - 0.36 

Q 01+U.8U 

u ' 3i -0.11 

n 90+0.U4 

U.28_ 0 .o6 

1.61 

178.4/161 


Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 ^ erg cm ^ s + All errors are in the 90% confidence level. 


due to the contamination of diffuse nonthermal X-rays 
around the source (in the SNR). We, therefore, investi¬ 
gate its 1.2- 7.0 keV spectrum with a single STEMS3D 
model, and obtain nH = 4.39 + 23 x 10 22 cm -2 , kT = 
0.51+02 keV, B = 9.73_ o . 59 xl0 14 G, A</> = 0.80++ 
/3 = 0. 25 + 03 , absorbed flux in 0.5-7.0 keV of 1.5 xlO -12 
erg cm -2 s _1 , and y 2 /dof = 121.4/98. Owing to the 
lowest flux and a large value of nH, we cannot obtain the 
upper error of B, which turns out to be the strongest 
surface magnetic field estimate among all sources. 

SGR 1900+14 is another source whose broadband 
spectra (1—200 keV) can b e mo deled with a BB plus 
a single PL component. IVrba et all (l2000lf reported that 
this source is embedded in a cluster of high-mass stars, at 
~ 12 — 15 kpc away. We fit the four XMM-Newton spec¬ 
tra of SGR 1900+14 in the 1.0—7.0 keV range in order to 
eliminate large deviations at the low energy most likely 
due to the absorption by the dense gas/dust in the clus¬ 
ter. We find that the surface temperature remains con¬ 
stant within errors among these observations (see Table 
ED- It is also not likely that the surface magnetic field 
strength varies over the course of about 2.5 yr of ob¬ 
servational span. We, therefore, link kT and B in our 
simultaneous fit so that our modeling would yield a com¬ 
mon value for these two surface parameters. We find that 
B = 8.0xl0 14 G, kT = 0.54 keV, /? ~ 0.11, and the twist 
angle is less constrained. If we further tie /3 and A <j> but 
allow the PL and STEMS3D normalizations to vary, we 
obtain the equivalently good fits (\ 2 /do/ = 361.9/392), 
with best-fit parameters of B = 7.81 xlO 14 G, kT = 0.56 
keV, A (f> = 1.92, and the (3 hit the lower limit (0.1) of pa¬ 


rameter space (Table El). These results suggest that the 
physical parameters do not evolve and SGR 1900+14 is 
a stable emitter. The ratio between the PL flux and the 
STEMS3D flux in 0.5-7.0 keV is - 30/70 (Figure!}. 

4.2.2. Variable Magnetars 

IE 1048.1-5937 and IE 2259+586 are two bright mag- 
netar sources, whose persistent X-ray flux varie d by a fac¬ 
tor o f 10 or more over t he timescale o f mon ths (Di b et al.l 
120091 : iZhu et all 120081) . iDib et all (1200911 performed a 
comprehensive timing analysis on ~ 10 yr RXTE data of 
IE 1048.1-5937, and showed that all three pulsed flares 
were accompanied by temporal events, e.g., significant 
pulse profile changes and/or a glitch. The infrared en- 
hance ment was also r eporte d at the onset of the 2001 
flare dWang fc Chakrabartvl l200lll . We fit all spectra of 
IE 1048.1-5937 with the single STEMS3D model (Fig¬ 
ure El) and find that the strength of the surface mag¬ 
netic field remains ~ 2.4xl0 14 G (see Group 1 in Table 
0}. By linking the B parameter for all observations, we 
find that the kT varies in the range of 0.39—0.54 keV, 
B = 2.42xl0 14 G, A <f> > 1.0, and /3 in the range of 
0.12—0.16. In order to provide constraints on the mag- 
netospheric parameters, in addition to B, we link the kT, 
A <j), and /3 for Obsl and Obs8, kT and /3 for Obs5 and 
Obs7, and kT only for Obs3 and Obs4. This grouping 
improves the fit quality slightly (x 2 /dof = 1113.9/915, 
Group 3 in Table Based on these results, we inves¬ 
tigated any possible correlations among the parameters 
of kT, /?, and 0.5-7.0 keV flux. We find that kT is cor¬ 
related with flux (the Spearman’s rank correlation co¬ 
efficient, p = 0.788, the probability of obtaining such a 







































4U 0142+61 


4U 0142+61 




Energy (keV) 


Energy (keV) 


RXS J170849.0-400910 



RXS J170849.0-400910 



Figure 6. Left panels: the X-ray count spectra of 4U 0142+61, 1RXS J170849.0-400910, and IE 1841-045 are fitted by the STEMS3D 
model ('Table |21) and the residuals are shown below. Right panels: the corresponding unfolded spectra. The dotted, dashed, and solid lines 
in the right panels mark the STEMS3D, PL components, and the sum, respectively. 
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Figure 7. Same as Figure|E]but for CXOU J171405.7-381031. 
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SGR 1900+14 SGR 1900+14 




Energy (keV) Energy (keV) 

Figure 8. Left panel: spectra of SGR 1900+14 are fitted with the STEMS3D model. Different colors correspond to the four different 
epochs of observations and their fitting residuals are shown in the bottom panel (Tabled}. Right panel: the corresponding unfolded spectra. 


Table 3 

Spectral parameters of SGR 1900+14 


Group 

Obs No. 

nH 

(10 22 cm- 2 ) 

kT 

(keV) 

o 

cq^ 

o 

A0 

(rad) 

P 

Flux 

X 2 /dof 


Groupl 

Obsl 

Obs2 

Obs3 

Obs4 

3 60 +U ib 

o.ou-oio 

U.48_o 04 

o.43i° 0 ;?l 

0-50 ±°;°i 
n 45+ 0 - 11 

io.o_ 0 6 
4-79 tUl 
8-95l°;^ 
8-99 2,3o 

U.ol_o.o8 
l-90_o .89 

1 qq+0.87 
l.uy_053 

2-0 1i4 

0 1 t)+ u,ua 

u - xo — 0.04 

0 20”*" 0 05 
U.ZU-003 

0 ifi+ 0 06 
u - 1D -0.04 

0 14+O.CM 

3.07 

3.08 

3.69 

2.93 

340.7/380 

Group2 

Obsl 

Obs2 

Obs3 

Obs4 

3 45 + u - i2 

0 54+ u - U2 

U - O4t -0.03 

8-OO^q js 

1.00_ 0> io 

i co+ 0.06 

J--t)0_oo5 

2 -°-o.i 

I 07+O.O8 
- 0.06 

0.10+ 002 

0.11+ 002 

0 12”*” 002 
u,lz -0.02 

0.10+ 002 

3.07 

3.09 

3.68 

2.93 

350.3/386 

Group3 

Obsl 

3 46~ l ~ u ' i;i 

0 5fi+ u - ui 

U. OO_0.04 

7.81^o.2i 

1 92 +u - u< ' 

0 io+o.° 4 

3.07 

361.9/392 


Obs2 






3.09 



Obs3 






3.68 



Obs4 






2.93 



Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 — ^ erg cm - ^ s — L All errors are in the 90% confidence level. 
Groupl: the only parameter of nH is linked for different observations in the simultaneous fit. Group2: parameters of nH, kT, 
and B are linked. Group3: parameters of nH, kT, B, A <j>, and /3 are linked. 
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correlation from a random data set, P = 0.020, and /? is 
anti-correlated with kT (p/P = —0.672/0.068), as well 
as the flux ( p/P = —0.783/0.021). It is important to 
note that all of these correlative behaviors emerge in all 
three steps of the fitting scheme presented in Table 0] 

The SGR-like bursts and a coincident timing glitch 
event occurred in the AXP IE 2259+586 on 2002 June 18 
(jKaspi et al. 120031) . The bursts and the glitch (that took 
place after Obs2 and before Obs3) were accompanied 
by the changes in the pulsed flux/profile and persistent 
flux/spectrum, providing the strongest ev idence of plas¬ 
tic deformation of the crust (|Woods et all 120041) . Hard- 
X-ray emission has bee n detected in IE 2259+586 by 
INTE GRAL and Nustar (|Kuiper et al.ll2006t lYoeel et ~ahl 
I2014D . Since the hard-X-ray tails contribute negligible 
flux below 7 keV, in which we are interested, their soft 
X-ray spectra are fitted with the single STEMS3D com¬ 
ponent (T ables0and Figure0. It had been suggested by 
IZhu et ahl (12008?) that the unabsorbed X'MM-Newton flux 
decay is best described by the PL and is correlated with 
the hardness. A simultaneous fit to all 10 observations 
yields a rather constant surface magnetic field strength. 
By linking B and refitting, we obtain B ~ 6.17xl0 14 
G (Group 2 in Table 0. To better constrain the 
magnetospheric parameters, we also merged some suc¬ 
cessive observations (namely, Obsl/Obs2, Obs4/Obs5, 
Obs7/Obs8, and Obs9/ObslO), which have similar sur¬ 
face emission properties (kT) and are either before or 
after the bursting episode. We find that the magneto- 
spheric twist angle rises to its maximum level follow¬ 
ing the bursting episode, and declines gradually over a 
timescale of two to three years. The particle velocity, 
on the other hand remains nearly constant around 0.2. 
Based on the results of Group 3 in Table 0 we find a 
strong correlation between A</> and flux with the Spear¬ 
man’s rank correlation coefficient of p/P = 0.942/0.005. 
This is an important indication of the twisted/untwisted 
magnetosphere scenario, which we will elaborate on in 
the Discussion section. 

4.2.3. Transient Magnetars 

XTE J1810-197 was discovered in 2003 when it 
suddenly brightened by about two orders o f mag¬ 
nitude above its quiesce nt flux (llbrahim et al.l 12004 
lHalpern fe Gotthelfl120051 ). As it is the first transient 
magnetar, this source has been visited with XMM- 
Newton more than 24 times and its spectra are we ll stud- 
ied, including investigations with STEMS (|G fiver et all 
I2007D . Due to the fact that the source shows no clear 
nonthermal component in quiescence, we only consid¬ 
ered the observations during the outburst decay before 
2005 March 18. There are two observations performed 
on 2003 September 8, and their spectral fitting results 
are consistent with each other. However, the S/N ra¬ 
tio of the second one (ObsID = 0161360401) is relatively 
low (due to its short exposure time of < 1 ks) and is 
excluded here. The surface temperature declines mono- 
tonically (from 0.52 ± 0.01 to 0.34 ± 0.01 keV) following 
the outburst of XTE J1810-197, while the magnetic field 
remains around (2.2 — 2.6) xlO 14 G, which is consistent 
with those found using STEMS modeling (|G fiver et al.1 
I2007D . In the meantime, the A (f> parameter typically re¬ 
mains around ~ 1.7 and /? reaches a lower limit of 0.1 
(Table 0. We also fit all four spectra with common val¬ 


ues of B and A <f>, which yields B = 2.52xlO 14 G and 
A <j) = 1.66. 

SGR 0501+456 was discovered because it exhibited 
a series of short energetic bursts de tected with multi¬ 
ple X-ray missions in 2008 August (jOdgiis et al.l l20lfl 
iKumar et al.ll2010h . The first XMM-Newton observation 
of this source took place while SGR 0501+456 was still 
burst active; therefore, there were a few hundreds of short 
bursts in about 45 ks of effective exposure. In order to 
exclude the times of these burst events, we constructed 
the source light curve in a 0.1 s bin size, and then used 
a cou nt rate cut-off criterion (< 55 cts s _1 ) dLin et all 
I2012H . SGR 0501+456 has also been detected in hard- 
X-ray band (lEnoto et al.ll2010a[) , but the hard-X-ray tail 
does not affect our fitting below 7 keV significantly. The 
soft-X-ray spectra are, therefore, fit with the STEMS3D 
model alone. We find that the surface magnetic field 
strength remains constant within errors (see Group 1 in 
Table 0. We refit all five spectra simultaneously to have 
a common value of B , and obtained it as B = 2.35xl0 14 
G. In this case, we obtain that the speed of magneto- 
spheric electrons declined from 0.2 to 0.15, while the 
twist angle showed a marginal evidence of variation be¬ 
tween 1.5 and 1.7 (Group 2 in Table 0. Among all of 
the evolving parameters, we find a significant correlation 
between only /3 and the flux, with p/P =1/0. 

Swift J1822.3-16 06 underwent an outburst in 2011 July 
(iLivingstone et aD 120111) . The timing analysis revealed 
a spin period of P = 8.4 s and a period derivative 
of P = 8.3 x 10 -14 s s _1 , which made this source 
the second magnetar with a_ low dipol e ma gnetic field, 
^timing — 2.7 x 10 13 G (|Rea et al.1 12012D . Its spec¬ 
tral behavior during the flux decay resembles those ob¬ 
served from XTE J1810-197 in the 2003 outburst, whose 
dipole magnetic field inf erred from spin parameters is 
^timing = 2.52xlO 14 G (IGamilo et al.ll2007lb The spec¬ 
tral evolution could be in terpreted as the thermal relax¬ 
ation in the NS crust dScholz et al.1120141) . The long¬ 
term monitoring data suggest that the temperature of 
the thermal component decreases slowly. Alternatively, 
the decreasing normalization, which is not presented in 
Table 0 indicates th at the emitting region shrinks durin g 
the outburst decay (IRea et all 12012b IScholz et al.ll2014D . 
In our investigation of this source with the STEMS3D 
model, we found that the surface temperature declined 
slightly from 0.54 ± 0.02 to 0.44 ± 0.02 keV in the first 
year of the outburst, and the surface magnetic field 
(B ~ 2.32xlO 14 G) and the electron velocity (/3 ~ 0.12) 
remained constant (Group 1 in Table 0 . Refitting all 
four spectra with a linked surface magnetic field and f3 
for all yields a significant fluctuation of A</> (see Group 2 
in Table 0 . We also fit all spectra with the same value 
of A (f> and the chi-square increases to 433.5/429, indicat¬ 
ing the a significant variation of A <j> with an F-test value 
lower than 0.005. 

IE 1547.0-5408 was discovered in 1980 
(lLamb fe Markertl 119811) , but has shown repeated 
X-ray flux enhancements in recent years: its flux went 
up by a factor of > 50 with the 2008 October outburst 
and further increased to levels of ~ 3 orders brighter 
than i ts lowest flux at the onset of its 2009 January out¬ 
burst dDib et al.l 120081 iKuiper et al.l 120121) . A transient 
hard-X-ray emission was also recorded with Suzaku 
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Table 4 

Spectral parameters of IE 1048.1-5937. 


Group Obs No. nH kT 

(10 22 cm -2 ) (keV) 


B 

(10 14 G) 


A<p 

(rad) 


/3 Flux x 2 /dof 


Group1 

Obsl 

1 19 +U.U 2 

- 0.01 

0 45 _l " u - u ' 5 

U.^O-o 03 

2 44 + u,i0 
Z< ^- 0.21 

1.90_i 21 

0 14+ U - U4 
u - x - 0.02 

0.56 

1060.6/902 


Obs2 


0 46 +0 ' 02 

9 qo + 0.08 

Z.OO -0 08 

l-96_ 0 os 

0 12+ 0 - 01 
U-1Z _ 0.01 

1.39 



Obs3 


0 42 +0,01 

U-4Z _ 0 Q1 

2 45 + 0 - 07 

1 4fi + 0,31 
4- ^°_o.09 

n 10 +O.O 2 
u - lo - 0.02 

0.96 



Obs4 


0 42+ 0,01 
u, ^ z - 0.01 

2-45l“;“ 

1 89 +0 ' 08 
x,oy —0.06 

0 16 - * -0 ' 01 
O-iO_ 0 . 0 i 

0.80 



Obs5 


0 40+ 0 - 01 
u.4u _ 0>01 

2 . 66 /g;“ 

d 

1 

O 

c4 

n 1 q+ 0-02 
U - 18 - 0.02 

0.72 



Obs 6 


0 53 + 0 - 01 

u.oo_o 01 

9 qf;+0.05 
z -'5O_0.04 

1 84 +0 ' 02 
1 -°^- 0.02 

0 i 2+ 0 - 01 

U-1Z — 0.01 

2.92 



Obs7 


0 40” 1 ” 0 ' 01 

9 cto+0-07 
z - oz ~0.05 

1 69” 1 ” 005 
i.oy _ 0 05 

0 16+ 0,01 
u ‘i°— 0.01 

0.54 



Obs 8 


0 4 Q+ 0-02 

U-4O 0 02 

9 qn+0.08 

z.ou 0 09 

1 QO+0.06 

t.uO-007 

0 lfi+ 0-02 
u - 1 D - 0.02 

1.01 


Group2 

Obsl 

1 1 Q + O.Ui 

1 ,lo - 0.01 

0 44 + u,U2 

u - 44 t - 0.02 

2 42+ U-Ui 

0.03 

-*■•■*■0-0.36 

n ifi+ u,ut> 

u - 1 D - 0.02 

0.56 

1108.3/909 


Obs2 


0 47 +O OI 
u -^'- 0 . 0 l 


1.96_o 05 

0 12 - *" 0-01 
U-l^_ 0.01 

1.39 



Obs3 


0 42“*" 0,01 
u.4z _ 0 01 


1 45+0-09 
1 , 4 O - 0.10 

0 14+0-01 
u-i^-o.oi 

0.96 



Obs4 


0 42+ 0,01 
u.4Z _ 0 01 


1 nn+0-07 

i.yu_oo7 

0 16+0- 01 
U-lO_ 0 . 0 l 

0.80 



Obs5 


0 3Q+0 - 01 

u.oy _ 0 01 


2 -°-o.i 

0 lfi+ 0-02 
u - 1 D - 0.02 

0.72 



Obs 6 


n 54 + 0,01 
u.o4 _ 0>01 


1 Rzl+ 0-02 

l-o4 _ 0 02 

0 12+ 0 - 01 
u-i^_o.oi 

2.91 



Obs7 


0 3Q+0 - 01 

u.oy _ 0 01 


1 54+0-07 

1-04—0 08 

n i 6+ 0 - 01 
u-io_ 0 .oi 

0.54 



Obs 8 


0 44+0 - 01 
u> ^- 0.01 


1 qi+ 0-07 
i - Ui -0.06 

0 16+ 002 
u - lo - 0.02 

1.01 


Group3 

Obsl 

1 1 Q“t“U-U^ 

■ LO — 0.02 

0 44 _I_U - Ui 
u *^- 0.01 

2 42+ u,U3 
-0.03 

1 06 +UiU 
1UD -o .10 

0 i 6 + u,u:i 
O-iO_ 0.01 

0.56 

1113.9/915 


Obs2 


0 47+0 - 01 
u -^'- 0.01 


1.96_ 0 os 

0 i 2+ 0 - 01 

U-1Z _ 0.01 

1.39 



Obs3 


0 42 - *" 0,01 
u.4z _ 0 01 


1 45+0 - 15 

1-40-004 

n 10 +O.OI 
u -lo_ 0 .oi 

0.96 



Obs4 


= obs3 


1 90 +0 08 

1.UU-0.64 

0 16 -1 " 001 
o-io _ 0 01 

0.80 



Obs5 


0 39 + 0 01 
u.o »_ 0 01 


to 

b 

l 

0 

0 16"*" 0 ' 01 
o-io _ 0 01 

0.72 



Obs 6 


0 ^4 +0 - 01 
u.o4 _ 0>01 


1 8 A+ 0-02 

l-o4 _ 0 02 

0 12+ 0-01 

2.91 



Obs7 


= 06 s 5 


1 45+0 - 04 
i-4O _ 0 12 

= o6s5 

0.54 



Obs 8 


= obsl 


= obs 1 

= obs 1 

1.01 


Note. — 

Flux: 0.5—7.0 keV absorbed flux 

n units of 10 — 

1 •] _9 _ 

erg cm s 

. All errors 

are in the 90% 

confidence level. Groupl: 


the only parameter of nH is linked for different observations in the simultaneous fit. Group2: both parameters of nH, and B are 
linked. Group3: /3 for Obsl and Obs 8 , kT and /3 for Obs5 and Obs7, and kT only for Obs3 and Obs4 are further linked. 


(lEnoto et all I2fll0bf ) and INTEGRAL (jKuiper et all 
I2012D during the decay of the 2009 outburst. The 
broadband X-ray spectrum can be described with a 
BB component plus a single PL component, and its 
spectrum evolved significantly in th e soft-X-ray s as well 
as in the hard-X-ray data ijKuioer et all 120121 . Here, 
we only use the XMM-Newton observation executed on 
2009 February 3, which is one week after the Suzaku 
observation (on 2009 January 28), and T of 1.41 is 
adopted based on th e 15-70 keV Suzaku HXD-PIN data 
(jEnoto et all l2010bl) . The low energy portion of its 
spectrum is contaminated because the s ou rce is located 
at the center of SNR G327.24-013 ( Gelfand fc Gaensler : 
120071 ). Therefore, we again fit the 1.2-7.0 keV spectrum 
in order to eliminate large deviations owing to the 
contamination of diffuse nonthermal X-rays at the low 
energy (Figure |TT]) . We find the best-fit parameters 
nH = 4.77l°;i 3 xlO 22 cm" 2 , kT = 0.50™ keV, 
B = 8.98™xl0 14 G, Acj) = 1.90™, (3 = 0 . 19 ™, 
absorbed flux in 0.5-7.0 keV F = 5.14 xlO -11 erg cm -2 
s -1 , and x 2 /dof = 96.8/102. Our modeling reveals that 
the STEMS3D component dominates the PL component 
below 7 keV (the flux ratio ~ 3 : 1). 

The transient SGR SGR 0418+5729 was discovered on 
2009 June 5 when its bursts triggered the Fermi/GBM, 
and then was reported as the fir st magn e tar w ith a low 
inferred dipol e magnetic fi e ld iR.ea et all (12010ft . Based 
on its spectra, IGiiver et al.1 (120111) argued that the value 


of the magnetic field at the surface is 10 14 G based on 
the STEMS fitting. The strong magnetic field case in 
this source was later confirmed by Tiengo et al. (2013). 
Nevertheless, the STEMS3D model cannot provide an 
acceptable fit to the spectrum of SGR 0418+5729 (ObsID 
= 0610000601). The canonical (BB+PL) model results 
indicate that the overall flux is dominated by a very hot 
thermal component (kT ~ 0.9 keV). The STEMS fitting 
results show that the particle velocity is /3 ~ 0.5 and the 
optical depth in the magnetosphere is r ~ 9. However, 
as can be seen in the Figure [3l the optical depth in the 
STEMS3D model is governed by both j3 and Act), and 
cannot be larger than two when the /3 is greater than 
0.3. 

The AXP CXOU J164710.2-455216 is lo cated in the 
youn g massive stars cluster, Westerlund 1 dMuno et al.1 
I2006D . The source underwent the first intense outburst 
on 2006 Septem b er and the se cond one on 2011 Septem¬ 
ber (iRodrfguez Castillo et al.ll2014D. Measu ring the spin 
period and period derivative, I An et all (12013!) suggested 
that the spin-down inferred magnetic field strength is less 
than 7 x 10 13 G. We attempted to fit six bright XMM- 
Newton pointing data during 2006 and 2011, and ob¬ 
tained statistically acceptable fits with y 2 /do/ ~ 0.99. 
However, the values of B in both STEMS and STEMS3D 
modeling peg at the low limit of 10 14 G; therefore, the 
obtained parameters are not reliable. We, therefore, do 
not include this source for further discussion. 

5. DISCUSSION 
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Table 5 

Spectral parameters of IE 2259+586. 


Group 

Obs No. 

nH 

kT 

B 

A (ft 

p 

Flux 

X 2 /dof 



(10 22 cm” 2 ) 

(keV) 

(10 14 G) 

(rad) 





Groupl 

Obsl 

0.82^;- 

0 39 +°- U3 
u,oy —0.05 

5 69 +u ' 88 

O.OU_o 43 

l- 34 -0.59 

fl 1 Q + 0.02 
U, 1 °- 0.02 

2.52 

871.1/928 


Obs2 


0 35 + 0 01 
u.oo-o oi 

6 ni +°- 26 
o,u± - 0.21 

1 45+0-04 

J-.4O-0 05 

0 20 + °' 01 

2.60 



Obs3 


0 33 +0 ' 01 
u.oo _ 0 oi 

9 61+0.30 
t7,u± —0.28 

i-97 _ 0 . 22 

0 22 +0 01 
0.01 

5.62 



Obs4 


0 36 +0 ' 08 
0.05 

6 14 + 0 . 9 4 
o,1 ^-1.50 

1 io + 0,46 
1U —0.26 

0 23 +0 ' 01 
u.z,o_o 01 

4.46 



Obs5 


0-42t°;°3 

5-79 ±811 

2 -°_0.4 

0 19+0- 02 
u,±y - 0.02 

4.59 



Obs6 


0-38l°;°? 

5.03l°;25 

1 90+ 0 07 
-L.yu-o 06 

0 20 + °' 01 
u.zu_o 01 

3.03 



Obs7 


n qo+0.03 
U.OO -0 02 

7 QO+ 1-01 

‘•* 2 - 0.75 

1 1 Q + 012 
±,±y -0.30 

0 20 + °' 01 
u-au_ 0 . 01 

2.96 



Obs8 


0 37 + 0 02 
u,0< - 0.02 

5 m +0 - 81 

o,Ui -0.13 

1-91-o.n 

0 20 + °' 01 
u.zu-o 01 

2.96 



Obs9 


0 37 + 0 - 02 
u - 0< -0.04 

5 38+ 0 - 21 

O.OO_0 53 

1-69-0.06 

0 20 + °' 01 
u-^u_o.oi 

2.94 



ObslO 


0 36 +0,04 
u - oo - 0.02 

5 87+ 0 - 51 
- 0.68 

1 46”*”°' 23 
1 -^°-0.24 

0 21 +0,01 
0.01 

2.87 


Group2 

Obsl 

0-79+"'"i 

n 41 + U - U2 

U.41_o 05 

6 17+ Uib 
°- 1 ' -0.13 

1 Q7 + u - 48 
1 * u ' -0.36 

fl 1 Q + 0.03 
U - 18 - 0.02 

2.52 

930.0/937 


Obs2 


0 37 + 0-01 
U.O/_o.oi 


1 oc+0.04 
J-ao_ 0 .ii 

0 20” 1 ” 0 01 
u - zu -0 01 

2.60 



Obs3 


0 4Q+0-01 

U. z tu_0.01 


2 -9_o. 1 

0 23 -1 " 0,01 
u.zo-o.oi 

5.63 



Obs4 


n qq+ 0-02 

U-o8_0.05 


1 in+0-30 
1 , 1 U - 0.12 

0 23“*" 0 01 
u.^o_ 0 .oi 

4.46 



Obs5 


n 43 + 0-02 

U.4:0-0 02 


2-0_o.3 

n 1 q+o .02 
U - 18 -0.02 

4.59 



Obs6 


0 37 + 0 01 
u -°' - 0.02 


1 30+0-09 

±.ou_oo9 

0 20” 1 ” 0 ' 01 
U.ZU _ 0 01 

3.03 



Obs7 


0 3fi+ 0 02 
U.OO_ 0>01 


1 45+0+3 
4:0-0.26 

0 20” 1 " 0,01 
u.zu_ 0 01 

2.96 



Obs8 


0 35 + 0 - 02 

u.oo-ooi 


1 50”*" 0,13 

±.ou_ 0 .26 

0 20” 1 ” 0 01 
uzu -0 01 

2.96 



Obs9 


0 37 + 0 01 
u, °' - 0.02 


1 24”*" 0,32 

J-a4_ 0>0g 

0 20” 1 " 0 01 
u.zu_ 0 01 

2.94 



ObslO 


0 37 + 0 01 
u,0/ - 0.02 


1 95+0-25 
±.zo_ 0 .27 

0 20” 1 ” 0 ' 01 
u.zu_o.oi 

2.87 


Group3 

Obsl/Obs2 

0-77++ 

Q 40+ U Ui 
U.4:U_0.01 

5 8 Q + U,U8 
°-oy_ 0.08 

1 4«+ u - U3 
J-4O_ 0 03 

n 1Q+ UUi 
u,iy -0.01 

2.56/2.60 

976.4/949 


Obs3 


0 43 + 0-01 
U.4:0_o.oi 


2-0_o.i 

0 21 +°- 01 
u.zi-o 01 

5.63 



Obs4/Obs5 


n 44 + 0-02 
u * 44 - 0.01 


1.96_ 0 13 

n 1 Q+o.oi 
b.lo_o.oi 

4.49/4.57 



Obs6 


0 41 +°-0i 


1 85”*"°' 07 
±.oo_ 0 .08 

0 17+0- 01 
u - 1 '-o.oi 

3.03 



Obs7/Obs8 


0 3Q + 0 - 02 

U.OU-o 02 


1 / 1 S+ 0.08 
1.4o_o.o4 

n iQ+001 
u,iy -0.01 

2.96/2.97 



Obs9/ObslO 


0 41 +°- 0 i 
u - 41 - 0.02 


1 oq-|-0.04 
J-ao_0.07 

0 20” 1 ” 0 ' 01 
u.^u_ 0>01 

2.94/2.87 



Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 — ^ erg cm - ^ s — + All errors are in the 90% confidence level. Groupl: the only 
parameter of nH is linked for different observations in the simultaneous fit. Group2: both parameters of nH, and B are linked. Group3: 
three pairs of successive observations (Obsl/Obs2, Obs7/Obs8, and Obs9/ObslO) are merged. 


Table 6 

Spectral fit results of XTE J1810-197 


Group 

Obs No. 

11 H 

( 10 22 cm” 2 ) 

kT 

(keV) 

B 

(10 14 G) 

A 4> 
(rad) 

p 

Flux 

X 2 /dof 


Groupl 

Obsl 

Obs2 

Obs3 

Obs4 

U.OO -002 

0 52+ u - Ui 
o.oz-o 01 

0 5 Q +0-01 

U.OU-002 

0 44+001 
u -44_o 01 

0 34+0 01 

U.O 4 _ 0 01 

9 99 +U.U 4 

0.04 

2 44 + 0-06 
z -^- 0.08 

2 50+0 06 

z.oU -006 

2 51+ 0-04 
zu± — 0.03 

1 66”*” U ' U2 
i.OO-002 

1.98_o.o5 

1 66”*"°' 09 
1 .OO -004 

1 54+0-03 
l-O4 _ 0 03 

O.IO + 001 

O.IO + 001 

0 11+0.01 
11 - 0.01 

0 1 o +0- 01 
U.1U_ 0 01 

3.76 

2.19 

1.32 

0.55 

559.7/442 

Group2 

Obsl 

Obs2 

0 Q3 + U,U2 

u.uo-ooi 

0-51++ 

0-49+°;°i 

9 C 9 + 0 .U^ 

^•O^_ 0 .02 

1 66”*” U ' U3 
i.OO-002 

O.IO + 001 

0 io+O.Q! 

3.74 

2.19 

750.7/448 


Obs3 

Obs4 


0 42”*” 0 01 

U.4Z _ 0 01 

0 32+ 0 01 
o.oz-o 01 



p p 

0 to 
+ 1 + 
000 
boo 

1.32 

0.55 



Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 — ^ erg cm - ^ s — + All errors are in the 90% confidence level. 
Groupl: the only parameter of nH is linked for different observations in the simultaneous fit. Group2: parameters of nH, B, 
and A <j> are linked. 
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IE 1048.1-5937 


IE 1048.1-5937 



Energy (keV) 


IE 2259+586 




Energy (keV) 


IE 2259+586 



Energy (keV) 

Figure 9. Same as Figure|6]but for IE 1048.1-5937 and IE 2259+586. 


Energy (keV) 


Table 7 

Spectral fit results of SGR 0501+456 


Group 

Obs No. 

nH 

kT 

B 

A <j> 

P 

Flux 

X 2 /dof 



(10 cm’^j 

(keV) 

(10 14 G) 

(rad) 





Groupl 

Obsl 

0-96+£oi 

0 44+ u,Ui 

9 90+e.uf 

z.zo—o 07 

i ,--,+0.06 

-i-.oi_o.o 8 

0 206"’" u,uut> 

U.ZUO_o 005 

3.72 

600.2/604 


Obs2 

0 40 + 0.01 
u. z ty_ 0 .oi 

9 qq + 0.07 
z.oO-007 

r q+0.04 

J-.oo_Q.08 

0 1 e;q + 0 - 008 

u.ioy_o 008 

3.13 



Obs3 


U.°U _ 0 .02 

9 34 + 0 - 06 
*-o4_ 0 10 

i '7i +0.05 
1 ‘ /1 -0.06 

0 1 50 +°.013 

U. 1 DU _ 0 021 

2.94 



Obs4 


0 50”*" 0 01 
U.OU—o 02 

2 41 +0.05 
z - 41 -0.08 

1 66”*" 0 ' 04 
l.OO—o 14 

0 14Q+ 0 - 010 

U.14O_0 016 

2.65 



Obs5 


0 44 - *" 0-01 
u> 44 - 0.02 

2 41 + 0 - 08 

1 65+ 0 08 
i.uo_o 30 

0 144”*" 0-013 
u - m _ 0.011 

1.44 


Group2 

Obsl 


0 45+ u - Ui 

u.4o _ 0 oi 

9 qc+O.Ui 

zo °- 0.04 

1 50 +u,uo 
i.ou_Q 06 

0 203"*” U ' UU4 

u.zuo _ 0 oio 

3.72 

605.9/608 


Obs2 


0 40 + 0.01 
u - 4 y _ 0.01 


i c-9+0.04 

i.oz_q os 

n i ko+o .006 
u.iao_o 006 

3.13 



Obs3 


0 50+ 0 01 

u.ou_o.oi 


1 71 +0-05 

A *' ± -0.06 

0 1 ^n+ 0.011 
U.lOU_o 019 

2.94 



Obs4 


0 49 + 0 - 01 
u -^ y - 0.01 


i fi5 + 0.04 

1 .00_o is 

0 149”*” 0 006 
u> 0.006 

2.65 



Obs5 


n 43 + 0 - 01 


1 65"*" 0,10 

i.oo_o 33 

n 140 + 0.012 
U.148_o 010 

1.44 



Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 — ^ erg cm - ^ s — ^. All errors are in the 90% confidence level. Groupl: 
the only parameter of nH is linked for different observations in the simultaneous fit. Group2: both parameters of nH, and B are 
linked. 
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0.5 1 2 5 

Energy (keV) 


SGR 0501+456 



Energy (keV) 



0.5 1 2 5 
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Figure 10. Same as Figure [6] but for transients XTE J1810-197, SGR 0501+456, and Swift J1822.3-1606. 


Table 8 

Spectral fit results of Swift J1822.3-1606 


Group 

Obs No. 

nH 

(10 22 cm" 2 ) 

kT 

(keV) 

B 

(10 14 G) 

A 0 

(rad) 

P 

Flux 

X 2 /dof 


Groupl 

Obsl 

Obs2 

Obs3 

Obs4 

0 42 +u - Ui 

u -^ z -0.01 

0 54 ~*~ u - uz 

U - O4l -0.02 

0 54+ 001 

U.04_q 02 

0 45 + 0 - 02 

U. 40 - 0.02 

0 44 + 0 - 02 
u -^ 4 -0.02 

9 *31 +O.U / 

ZlO1 -0.08 

Q qq + 0.06 

Z.OO- 0.04 

2 34 + 0 - 08 

z -°4_o .08 

9 qcc+ 0.08 

z.oo-o 09 

1 49 +U.U 4 

J-- 4 Z- 0.07 

1 C9+0.04 

^•O^_o.o6 

J--oU_o.o3 

1 cn+ 0.05 
i.oy 0 04 

0.12+ 0 ' 03 

O.II + 003 
n 10+O.OI 
u,lc> -0.01 

0 12+ 001 
u - lz -0 01 

1.90 

1.31 

0.43 

0.28 

415.5/420 

Group2 

Obsl 

Obs2 

Obs3 

Obs4 

0 42 +u,Ui 
u -^ z -0.0i 

0 54 + u - Ui 

u.o 4 _ 0 oi 

0 153+ 001 
u.oo-o.oi 

0 4 fi+ 0 - 01 
u - 40 _o 01 

0 43+ 0 01 
u - 4 o_o. 01 

9 q 9 + U.U4 

1 42+ u - U4 
— 0.07 

1 ci +0.04 
0.06 

1 oq+ 0.04 

t-oU_o.o3 

1 cq+ 0.05 
i.°y _ 0 .04 

0 12 +U,ui 
u - iZ -o.oi 

1.91 

1.31 

0.44 

0.28 

420.5/426 


Note. — Flux: 0.5—7.0 keV absorbed flux in units of 10 — ^ erg cm - ^ s — b All errors are in the 90% confidence level. 
Groupl: the only parameter of nH is linked for different observations in the simultaneous fit. Group2: parameters of nH, B, 
and /3 are linked. 
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Figure 11. Same as Figure[6]but for the observation of IE 1547.0-5408 taken on 2009 February 3. 
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Figure 12. Upper panel: nHgTBMS3D vs. nH( BB+PL ). Bottom 
panel: -Bg TEMS3D vs. -Bstems- 


Table 9 

Log of the Magnetic Field 


Source 

-^timing 

(10 14 G) 

-BstEMS 
(10 14 G) 

®STEMS3D 

(10 14 G) 



4U 0142+61 

1.3“ 

5 45+ u - 2b 

°- 40 — 0.14 

5 «7+ u - bS 

-0.49 

1RXS J170849.0-400910 

4.6 b 

5 QQ + 0 - 32 

o.yy_o 56 

6.03l“' 3 ° 

IE 1841-045 

6.9 b 

3A2f° 0 f 0 

4-i6±g;lS 

CXOU J171405.7-381031 

5.0 C 

6-99t°j 4 

9-73 _ 0 . 59 

SGR 1900+14 

7.0^ 

4 Q5+0.64 

4.y°_ 0< 22 

1 oi +0.61 
' - 81 - 0.21 

IE 1048.1-5937 

3.9 e 

o 97 +O.O 4 
-0.04 

2 49+ 003 
0.03 

IE 2259+586 

0.59 f 

5 02~*” 0 ' 13 
o.uz — o 07 

5 sq +°- 08 

o.oy_ 0- 08 

XTE J1810-197 

2.19 

2 48 +0 ‘ 02 
z -^°- 0.02 

9 ^9+0.62 

SGR 0501+456 

1.9 h 

9 -i 9+0.04 
^•J-^-0.04 

z -°°_0.04 

Swift J1822.3-1606 

0.135* 

9 1 q+0-06 

^•- lo —0.04 

9 09 +O.O 4 
04 

IE 1547.0-5408 

3.29 

7.95-2.54 

8.98/g; 3 6 


Note. Period derivatives of magnetars are highly 
variable, and the derived Btiming could be significantly dif¬ 
ferent from the values listed here. 


(a - ) IDib et al.l 120071: fbl IDib et all 120081: (cl ISato et all 
l2Q10l; (df IMere ehetti e t al.1 1200a fel IP ib et a l.1 120091 : 


ff) Gavriil & Kaspi :2004" (g) Camilo et abj [2007 

; (hi 

Goghs et al. 20l(| (i) Scholz et al.i2014j; (j) Dib et al. 

2012 . 


In this paper, we developed a physically motivated nu¬ 
merical model to account for the magnetar surface emis¬ 
sion and magnetospheric scattering in three dimensions. 
We applied our model to the X-ray spectra of 13 magne¬ 
tars, including persistent X-ray emitters as well as tran¬ 
sient sources. We find that most spectra analyzed here 
were adequately described with the STEMS3D model, 
yielding important physical parameters, such as the sur¬ 
face magnetic field strengths and degree of twisting in 
the magnetosphere. To compare with our results, we also 
modeled all magnetar spectra with the STEMS model, 
which employs the same atmospheric emission proper¬ 
ties, but incorporates an only ID s catter ing scheme on 
the emergent radiation (iGiiver et al.ll2007T ). We find that 
the surface magnetic field strengths obtained by fitting 
the STEMS3D model are generally slightly higher, but 
are still in agreement within errors with those obtained 
with the STEMS model (see the lower panel of Figure 

m i- 

Dipole magnetic fields inferred from spin parameters 
(^timing) in 7 out of the 11 sources are consistent with 
the values obtained from the STEMS3D model fitting 
(Table ©. The STEMS3D model does not provide 
an acceptable fit to the spectrum of the first low-P 
magnetar, SGR 0418+5729, for reasons as already dis¬ 
cussed in §4.2.3. While for another two magnetars with 
^timing < 10 14 G, IE 2259+586 and Swift J1822.3-1606, 
the strengths of B given by the STEMS3D model are 
about one order of magnitude higher than their dipole 
magneti c fields. However, taking both twisted mag neto- 
spheres ([Thompson et al.l l2d02tlBeloborodovl 120091 ) and 
particle outflows ( Tong et al.ll2013T) into account, the ex¬ 


pected spin-down torque (P) will be larger than a pure 
dipole case with the same poloidal magnetic field. The 
conflict between larger values of PstemS 3 D and timing 
results could be reconcil ed by e ith er a small magnetic 
inclination angle (e.g.. lTong &; Xull2012l) and/or a more 
complex structure of magnetosphere (i.e., higher order 
multipoles). In a multipolar field, the strong surface 
magnetic field decays more rapidly with r adius; thus, a 
too fast spin-down rate can be avoided (jTurolla et alJ 

1ml). 

On the other hand, an important caveat for the in¬ 
ferred dipole field is that, magnetars are very noisy and 
their period derivativ es can vary sign i ficant ly within a 
short timescale (e.g., lArchibald et al.l f2015l and refer¬ 
ences therein,), leading to fluc tuations in -Burning , which 
is unphysical. For example, iLivingstone et all ([2011D 
measured a spin-down rate of P = 2.54 x ICE 13 s 
s” 1 using data collected in the first four months fol¬ 
lowing the 2011 July outburst of Swift J1822.3-1606. 
iRea et ahl d2013f) estimated P = 8.32 x 10 -14 s s -1 with 
the data covering 2011 July - 2012 April, and recently 
iScliolz et al.l (120141) took the new Swift data into account 
and reported P = 2.15xl0~ 14 ss _1 . As a consequence, 
the corresponding Ptiming decreased from the previous 
value of 4.7xl0 13 G to 1.4xl0 13 G listed in Tabled 
The decreasing spin-down r ate can be interpreted as ei¬ 
ther a sig n of wind braking dTong fc X u 2013) or a glitch 
recovery dScholz et al.ll2014l) . Our spectral investigations 
with the STEMS3D, however, yields magnetar-like sur¬ 
face magnetic field strength (2.3xl0 14 G), which is likely 
the case given that, at the very least, the source exhibited 
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energetic bursts. 

We find that the temperatures of the magnetar surfaces 
vary between 0.24 and 0.56 keV, and the magnetospheric 
electrons are non-relativistic (/? < 0.2). We investigated 
multiple XMM-Newton observations typically spanning 
years of five magnetars and find that the strength of the 
surface magnetic fields remains nearly constant. We also 
modeled all magnetar spectra with the phenomenologi¬ 
cal (BB+PL) model to better establish the interstellar 
absorbtion aspects toward magnetars studied here. It 
is well known that the simple PL model increases very 
rapidly at low energies, and as a consequence, a larger 
absorption column density is required to account for the 
overestimated PL flux, especially for the steep PL indices 
llDurant fc van Kerkwiikn2006t iGiiver et al.ll2008h . In all 
sources except IE 1841-045, the hydrogen column density 
obtained from the STEMS3D model fits (nHsTEMS 3 D) is 
lower than the value derived from the (BB+PL) model 
fits (nH (BB+ p L) ; Figure fl2l). We argue that, given the 
empirical nature of the (BB+PL) model, the density of 
interstellar hydrogen obtained with the STEMS3D is a 
more reliable indicator. 

As already stated, the surface magnetic field values 
from the STEMS3D model mostly agree with those ob¬ 
tained with the STEMS model results (Table 0 . How¬ 
ever, the STEMS model makes simplifying assumptions 
about the magnetic field geometry, i.e., treating the scat¬ 
tering region as a plane-parallel slab; therefore, it cannot 
offer information about the geometry of the magneto¬ 
sphere. Alternatively, the 3D self-similar magnetic con¬ 
figuration is characterized by a free parameter, A</, in 
the STEMS3D model. Our fitting results indicate that 
the magnetospheres in most magnetars are highly twisted 
(A</> >1). It is worth noting that we also included a flat 
PL (P = 1.32) component in the fitting for IE 1841- 
045, which could be the reason why it has a large value 
of n7?sTEMS3D- The derived small magnetospheric twist 
angle Acf> = 0.31 +p.ii and large /3 = 0.28 +q;q B should be 
treated with caution, and longer exposure observations 
are required to check the magnetospheric parameters of 
IE 1841-045 in the future. 

Compared with the STEMS model, the STEMS3D 
model embodies some other significant improvements: 
(1) in the ID geometry, the magnetic field is fixed to 
follow a B oc r~ 3 dependence, and the optical depth r is 
independent of the electron velocity /?. In contrast, the 
radial dependence of the magnetic field is consistently ob¬ 
tained from A</, while the angle-averaged optical depth 
r is determined by /3 and Aft. (2) The charged particles 
are assumed to have a top-hat velocity distribution in the 
ID model, while a more realistic distribution, involving 
the effects of bulk motion and also thermal velocity dis¬ 
tribution of electrons, i s used i n th e STE MS3D model. 
Following the work in iNobili et al.l (120081) , we adopt the 
ID Maxwellian distribution at the electron temperature 
T e , which is further linked with /3 by assuming equipar- 
tition between thermal and bulk kinetic energy in order 
to reduce the number of parameters in the model. (3) 
In the 3D configuration, the optical depth varies signifi¬ 
cantly with angles and vanishes at magnetic poles in the 
twisted magnetosphere. Thus, it is possible to detect the 
surface radiation which is mostly carried by extraordi¬ 
nary mode photons, if one of the magnetic poles is in our 


line of sight. The polarization and the proton cyclotron 
absorption lines a re expected to vary with the rotational 
phase. Note that iTiengo et all (I2013D recently reported 
a phase-variable absorption feature observed during the 
outburst of SGR 0418+5729. 

We also aimed to understand the nature of variable 
and transient magnetars by investigating their spectral 
evolution. For this purpose, we analyzed multiple XMM- 
Newton observations of six magnetars. IE 1048.1-5937 
has been quite active in the past 10 yr: three X-ray 
flares w ere recorde d with RXTE in 2001, 2002, and 2007 
March dDib et all 12009 1). The eight XMM-Newton ob¬ 
servations we used here are sampled before, after and 
in between these flares: Obsl was performed before the 
first flare; Obs2-5 were executed between the second and 
third flares; Obs6-7 were performed after the third flare; 
comparing with Obs7, the source brightened by a fac¬ 
tor of two in Obs8 (on 2013 July 22, after the end of 
RXTE mission). We find that the surface temperature 
monotonically decreases with declining flux from Obs2 to 
Obs5 and from Obs6 to Obs7 (see Group 1-3 in Table 0). 
This is likely due to cooling of the NS crust that is heated 
by the mechanism leading to these flares. However, the 
magnetospheric twist angle cannot be constrained in this 
source. 

In the magnetar f ra mewor k, both the crustal cool¬ 
ing (e.g., iPons fc Rea I l2012f) and the magnetospheric 
relaxation ie.g.. iBernardini et, al.ll2009fl models are pro¬ 
posed for the flux relaxation in different sources. The 
dense XMM-Newton observations on IE 2259+586 pro¬ 
vide a good opportunity to study this topic. The first 
two data sets are observed before the burst event, and 
the other eight observations monitored the flux relax¬ 
ation of the outburst in the following three years. The 
STEMS3D fitting results suggest a large magnetic field 
~ 5.9xlO 14 G for the source, and the A(j> suddenly in¬ 
creased to a maximum from Obs2 to Obs3 then decreased 
following the flux decays. The only exception, Obs4, 
which was executed ~ 20 days after the burst, has a 
small A(f> = 1.1. Note that, Obs4 and Obs5 were per¬ 
formed on the same day; however, the latter yields a 
large A <j> ~ 2. Such a small discrepancy between Obs4 
and Obs5 was a lso seen in the (BB+PL) fitting results 
(IZhu et al.ll2008l b i.e., the flux ratio between PL and BB 
components varied from 2.6 (Obs4) to 2.1 (Obs5). It 
is not likely that the magnetosphere could untwist and 
retwist in a few days, and the small A <j> derived from 
Obs4 is dubious. Group 2 in Tabic [5] shows a strong 
correlation between A<p and flux with the Spearman’s 
rank correlation coefficient of p/P = 0.661/0.038. Since 
the Spearman’s rank correlation coefficient does not con¬ 
sider the parameter uncertainties, it would underesti¬ 
mate the significance level due to the spectral fluctua¬ 
tions. Therefore, we merge three pairs of observations 
(Obsl/Obs2, Obs7/Obs8, and Obs9/ObslO) that have 
similar spectral parameters and flux levels, as well as 
Obs4/Obs5, and the correlation between Acj) and flux 
becomes more significant ( p/P = 0.942/0.005, Group 3 
in Tabled . When A <f> is linked, the x' 2 increases by more 
than 51.9, which suggests the significant variation of A <f> 
with an F-test probability lower than 2 x 10^ 9 . There¬ 
fore, our finding of an increased twist angle accompanied 
with the outburst (occurred between Obs2 and Obs3) 
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and then a decrease following the source decays, agrees 
perfectly with t he twiste d/untwisted magnetosphere sce¬ 
nario ([Thompson et al.ll2002ll . 

The STEMS3D fitting results of both XTE J1810-197 
and Swift J1822.3-1606 show that the surface temper¬ 
ature monotonically declines with decreasing luminos¬ 
ity. The former source ha d been studied in detail with 
the STEMS model (IGiiver et al.l I2007h . On the other 
hand, it is the first time that we evaluated the sur¬ 
face magnetic field of Swift J1822.3-1606 via continuum- 
fitting, and obtained a significantly stronger magnetic 
field than the dipole magneti c fi eld inferred from the 
spin-down rate (|Scholz et al.1 I2014D . It is interesting 
to note that XTE J1810-197 and Swift J1822.3-1606 
share t he similar spectral sequences followi ng their out¬ 
bursts (|R,ea et al.ll2012t IScholz et, al.ll2014H . Neverthe¬ 
less, the timing analyses infer a typical magnetar mag¬ 
netic field for XTE J1810-197 (B timing ~ 2.1 xlO 14 G, 
iCamilo et al.ll2007h while a low value for the latter source 
(-Btiming ~ 1.4xl0 13 G, IScholz et al.ll20l4) . Our spec¬ 
troscopic measurements suggest a similar magnetic field 
strength of -BstemS3D ~ 2.3xl0 14 G for these two tran¬ 
sients. The monotonic decline of surface temperatures 
during the outburst decays in both XTE J1810-197 and 
Swift J1 822.3-1606 are remarkable evidence o f crustal 
cooling (IGiiver et al.1120071 : IScholz et al.1I2014D . Addi¬ 
tionally, we also find a significant variation of A <f> in 
Swift J1822.3-1606 (with a confidence of 99.5%). We 
also find a decline of the magnetospheric electron ve¬ 
locity 0 in the first ~ 40 days of the 2008 outburst of 
SGR 0501+456 (Table [TJ) , implying the energy dissipa¬ 
tion of charged particles, probably by the radiative drag. 

Different scenarios have been proposed to ac¬ 
count for the har d-X-ray e missio ns from magnetars. 
iThompson fc Beloborodov! (I2005f) put forward two 
mechanisms that soft-q-rays either arise from the ther¬ 
mal bremsstrahlung emissions in a thin surface layer 
heated by a returning current or t he syn chr otro n emis¬ 
sions from pairs (see also iBeloborodov fc Thompson! 
I2007D . Alternatively, iBaring fc Hardind ( 20071 ) sug- 
gested that the high energy tail was produced by the 
resonant up-scattering in the magnetosphere. The 3D 
Monte Carlo simulations show that the soft-X-ray pho¬ 
tons can be up-scattered to ~ 100 keV if the scattering 
particles are energetic, i.e., 0 > 0.6 and the corresp ond¬ 
ing e l ectron temperatu re kT e > 100 keV (N obili et al.l 
120081 : [Zane et alJI2009h . However, our spectral analysis 
shows that the A (f> > 1 and the 0 ~ 0.2 in most cases, 
while flatter spectra and stronger pulsations detected at 
the hard-X -ray clear ly point to a different population of 
electrons. iBeloborodovI (120131 ) proposed that e 41 pairs 
discharge near the neutron star and create the relativis¬ 
tic outflow that further scatters thermal photons to high 
energies. This model provided a satisfactory explanation 
for the hard- X-ray luminosity, spectral slopes, and p ulsed 
profiles (e.g.. lHascoet, et afll2014t iVogel et, alJl2011l) . and 
it may also contribute an non-negligible flux in the soft- 
X-ray band. 

It should be noted that the STEMS3D model still 
has some limitation s, as other magne tosp heric scatter¬ 
ing m odels do (see I Zane et all 120091 and Bcl oborodovl 
l2013f l. Both spectral and temporal analyses have demon¬ 
strated that the thermal emission is not uniformly dis¬ 


tributed on the surface of magnetars but are confined 
in a small region (e.g., lOzell 120(5 lAlbano et al.l l2010h . 
However, because the uniform distribution is adopted 
for the seed photons, the STEMS3D model cannot 
provide an acceptable fit to the spectra of magnetars 
in quiescence, which lacks nonthermal emission. As 
the transients return to the quiescent state, the spec¬ 
tra of the three transients studied in our paper (i.e., 
XTE J1810-197, SGR 0501+456, and Swift J1822.3- 
1606) became significantly softer than those at outburst 
levels, and the second or even t hird thermally emitting 
surface regions were reported dBernardini et al.l 120091 : 
iGamero et al.ll201~it IScholz et a. 1 .1120141) . According to the 
(BB+BB) or (BB+PL) fits, it is found that the thermal 
emitting area decreases following the outburst decay. In 
these cases, the distinctive spectra would be formed for 
the following reasons: (i) In the 3D twisted magneto¬ 
sphere, the optical depth and upscattering efficiency vary 
with latitude because the current density and direction 
(along the magnetic field line) depend on the position. 
The optical depth has a maximum and a minimum at 
the equator and the poles (no current exists along the 
magnetic axis), respectively. Thus, a hot spot presented 
at the pole would gain less energy from particles in the 
magnetosphere, resulting in a soft spectrum, (ii) Pho¬ 
tons would be scattered away from the initial propaga¬ 
tion directions and might move into or out of the sight of 
view. The viewing effect should be taken into account, 
(iii) The superposition of two thermal components at 
different temperatures modifies the spectral profile fur¬ 
ther. To date, however, it is impossible to incorporate 
all of these effects since it will involve too many parame¬ 
ters. We emphasize that those extremely soft spectra at 
very low flux levels cannot be fitted by the STEMS3D 
model nor the other normal Compton scattering models 
(e.g., ID RCS, 3D RC S, STEMS models). Neverthe¬ 
less, I Zane et all (1200911 claimed that a softer spectrum 
is caused by a non-homogeneous distribution of sur face 
emissi on (effects (i) and (ii)). Meanwhile. iTurolla et alj 
(201 I i) also suggested that a surface thermal component 
dominated spectrum can be reproduced by considering 
eff ects (i) and ( iii). 

IBeloborodovI (120091) argued that the puzzling behav¬ 
ior of XTE J1810-197 can be explained with a un¬ 
twisting magnetosphere, which is divided into a current- 
free (“cavity”) region and a current-carrying bundle (“j- 
bundle”) of field lines. The hot spot in the polar re¬ 
gion is interpreted as the footprint of the j-bundle, which 
shrinks with time. As the magnetosphere untwists, the 
cavity expands and the j-bundle becomes so narrow that 
only a small fraction of the surface radiation can be scat¬ 
tered, producing very soft spectra in quiescence. How¬ 
ever, in this model it is difficult to explain that lack of 
flux enhancement whe n glitches occurred i n the stable 
sources. Alternatively, iPons fc Rea I (|2012l ) interpreted 
transient behaviors as the magneto-thermal evolution, 
and suggested that short X-ray bursts and glitches might 
always be accompanied by a flux variation; neverthe¬ 
less, the long-term flux of bright/stable magnetars (e.g., 
4U 0142+61) cannot be enhanced significantly. 

In some cases, magnetic configurations might be more 
complex than a dipole, and multipoles fields are re¬ 
quired to avoid too fast spin-down as we discussed above. 
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In addition, the theoretical calculation indicates that 
the output luminosity from a globally twisted magne¬ 
tosphere with A (f> > 1 would be at the or der of 10 36 erg 
s~\ which is higher than the observed dThomoson et, all 
120021 : lBeloborodovll200'9li . There are now some observa¬ 


tional phenomena in favor of partially, i nstead of glob¬ 
ally, twisted mag netospheres (e.g., iWoods et all l2007t 
iPavan etTa l. 2009), w hich ha ve been considered in the lit- 
erature (iBeloborodovI 120 09). Finally, although we make 
an assumption of constant velocity for charged particles, 
the particles should be decelerated via scattering of pho¬ 
tons, broadening their velocity distribution. 
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